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ABSTRACT 

Accurate photometric reconstruction of the Lunar surface is im- 
portant in the context of upcoming NASA robotic missions to the 
Moon and in giving a more accurate understanding of the Lu- 
nar soil composition. This paper describes a novel approach for 
joint estimation of Lunar albedo, camera exposure time, and pho- 
tometric parameters that utilizes an accurate Lunar-Lambertian re- 
flectance model and previously derived Lunar topography of the 
area visualized during the Apollo missions. The method intro- 
duced here is used in creating the largest Lunar albedo map (16% 
of the Lunar surface) at the resolution of 10 meters/pixel. 


1. INTRODUCTION 



High resolution and accurate albedo maps of planetary surfaces in 
general, and the Lunar surface in particular, are used in landing site 
selection, mission planing, understanding soil composition, and as 
educational resources. Accurate albedo maps effectively support 
landing operations as they are used to estimate the landing site 
appearance under the illumination conditions (Sun and spacecraft 
position) at landing time. In planetary science the albedo maps 
advance our understanding of soil porosity, surface composition 
and formation through correlations with terrain models and other 
data mapping (multi-spectral, mineralogical) products. 

This paper describes a method for albedo map reconstruction 
using imagery retrieved by the Apollo 15, 16 and 17 missions. The 
Apollo Metric Camera was flown on an orbit at approximately 100 
km above the Lunar surface. It was a wide field of view (75°) or- 
bital mapping camera that photographed consecutive images with 
approximately 80% overlap. The digital scans of these film im- 
ages [9, 3] capture the full dynamic range and resolution of the 
original film, resulting in digital images of size 22,000 x 22,000 
pixels, representing a resolution of 10 meters/pixel and covering 
approximately 16% of the Lunar surface. Each pixel of the Apollo 
Metric Camera images was formed by a combination of many fac- 
tors, including albedo, terrain slope, exposure time, shadowing, 
and viewing and illumination angles. The goal of albedo recon- 
struction is to separate the contributions of these factors. This is 
possible in part because of redundancy in the data; specifically, 
the same surface location is often observed in multiple overlap- 
ping images. The albedo reconstruction technique includes all of 
the above quantities in an image formation model and re-estimates 
them from the existing imagery. The overall system is shown in 
Figure 1 . The components of this block diagram are explained in 
more detail in the following sections. 


Fig. 1. The overall albedo reconstruction system. 


2. PHOTOMETRIC MODEL 

Starting with the first images from the Apollo missions, a large 
number of Lunar reflectance models were studied [6, 7, 8]. In 
this work the reflectance is derived from the Lunar-Lambertian 
model [6, 5]. As shown in Figure 2, we define the following unit 
vectors: n is the local surface normal, and 1 and v are directed 
at the locations of the Sun and the spacecraft, respectively, at the 
time when the image was captured. We further define the angles i 
separating n from 1, e separating n from v, and the phase angle a 
separating 1 from v. 

The reflectance model used in our approach is given by 
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where L(a) is a weighting factor between the Lunar and Lamber- 
tian reflectance models [2] that depends on the phase angle and 
surface properties. R is a photometric function that depends on 
the angles a, i and e. The parameters c\ and C 2 play a significant 
role in modeling the surface reflectance at low phase angles. This 
can be noticed in the reduction of the banding artifacts in Figure 8 
on the top right. 











Fig. 2. Illumination and viewing angles used by the Lunar- 
Lambertian reflectance model. 


3. PARAMETER ESTIMATION 

The computation of the reflectance (”R Map” in Figure 1) using 
Equation 1 is central to our albedo reconstruction technique. Let 
Iij , Aij, Rij be the Apollo Metric Camera image value, albedo 
and reflectance at pixel ij , and T be a variable proportional to the 
exposure time of the image. Then 

Iij — T Aij Rij. (2) 

Note that the image formation model described in Equation 2 does 
not take into consideration the camera transfer function since the 
influence of the non-linearities of the camera transfer function plays 
a secondary role in the image formation model [2]. From Equa- 
tion 2 it can be seen that when the observed pixel value, expo- 
sure time, and reflectance value are known, the image formation 
model in Equation 2 provides a unique albedo value. However, 
these observed values are subject to errors arising from measure- 
ment (exposure time), image scanning process (image value) and 
stereo modeling errors (reflectance), resulting in imprecise albedo 
calculations. The method presented here mitigates these errors by 
reconstructing the albedo of the Lunar surface from all the over- 
lapping images, along with their corresponding exposure times and 
DTM information. The albedo reconstruction is formulated as the 
least squares problem that minimizes the following cost function 
Q: 


Q = J2J2 [Aj - (3) 

k ij 

where the superscript k denotes the variables associated with the 
k - th image and Sfj is a shadow binary variable, Sjfj = 1 when the 
pixel is in shadow and 0 otherwise. The weights w\j are chosen 
such that for each image k they have linearly decreasing values 
from the center of the image (w- = 1) to the image boundaries 
(■ Wij =* 0). The choice of these weights insures that the recon- 
structed albedo mosaic is seamless. The S'- values play an impor- 
tant role in accurate albedo estimation by discarding image pixels 
that are in shadow and for which the DTM and the reflectance 
models are unreliable [1, 4]. The DTM used in computing the 
reflectance value R . - at each pixel is obtained by a weighted aver- 
age (using the above weighting coefficients) of the DTM generated 
using stereo from all pairs of consecutive Apollo Metric Camera 
images. Figure 3 shows the DTM before and after the weighted 
DTM averaging and denoising process. It can be seen that the 


noise artifacts in the original DTM are reduced in the averaged 
DTM while the edges of the large craters and mountain regions are 
well-preserved. The color shaded averaged DTM of the Apollo 15, 



Fig. 3. Hill- shaded maps generated using (left) single local DTM 
and (right) denoised DTM derived from weighted average of over- 
lapping DTMs. Our denoising approach preserves the structure 
while reducing the artifacts shown in the insets. 

16 and 17 missions is shown in Figure 4 over the image mosaic of 
the Clementine mission. 



Fig. 4. Color shade of the enhanced digital terrain model of the 
Apollo zone. 

The optimal albedo reconstruction [10] from multiview im- 
ages and their corresponding DTM is formulated as a minimiza- 
tion problem of finding 

{Aij , f k , ci } = arg min Q (4) 

Aij,Tk, Cl 

for all pixels ij and images k , where Q is the cost function in 
Equation 3. An iterative solution to the above least square problem 
is given by the Gauss-Newton updates given below. 

• Step 1 (initialization): Compute the enhanced DTM as de- 
scribed earlier and the weights. Normalize the weights so 
that the sum of weights over all images is 1 at each pixel. 
Initialize the exposure time as inversely proportional to the 
average image reflectance. Initialize the phase coefficients 
ci and C 2 to some reasonable values. Initialize the albedo as 
the arg min of the cost function Q for fixed exposure time 
and phase coefficients 
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• Step 2: Re-estimate the exposure time using 
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• Step 3: Re-estimate the phase coefficients using 
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and analogously for C 2 . The partial derivatives dR^j / dc\ 
and dRij / dc 2 are computed from Equation 1 . 

• Step 4: Re-estimate the albedo using 
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ring around its rim in the reconstructed albedo image (Figure 7 
top-right). Figure 7 (bottom) shows the illumination removal effect 
in the reconstructed albedo for the Apollo 15 landing site. The 
illumination/shadow removal effect as well as the seam reduction 
is shown in Figure 8 (compare the left and right images). Pivotal 
to reducing the seam artifacts in the top-right image in Figure 8 
was the use of the phase coefficients c\ and C 2 in the reflectance 
formula, Equation 1 . 

The reconstruction error map (Figure 6) is computed at each 
pixel using Equation 9. Brighter values have higher estimated er- 
rors, which in turn implies a larger discrepancy between the pixel 
values at that location in the original overlapping images. It can 
be noticed that larger albedo reconstruction errors occur around 
crater ridges and are due to larger terrain reconstruction errors in 
these regions. Larger than average albedo reconstruction errors 
are also noticed at boundaries between orbits on a East- West Lu- 
nar traverse and are due to the small fraction of overlap between 
images taken in adjacent orbits compared to the overlap between 
consecutive images within the same orbit. 


Step 5: Compute the error cost function Q for the re-estimated 
values of the albedo and exposure time. 

Convergence: If the error between consecutive iterations 
falls below a fixed threshold, then stop. Otherwise return to 
step 2. 

Step 6: Estimate the error map of the albedo reconstruction 
at each pixel using the formula 

- - 

4. EXPERIMENTAL RESULTS 



The main results of this research are the reconstruction of the Lu- 
nar albedo from images taken by the Apollo Metric Camera flown 
onboard the Apollo 15, 16 and 17 missions, and the correspond- 
ing error map. The reconstructed albedo (Figure 5) covers ap- 
proximately 16% of the Lunar surface at an unprecedented resolu- 
tion of 10 meters/pixel for this coverage. For comparison, LRO- 
WAC albedo mosaic covers 100% of the Lunar surface at 100 me- 
ters/pixel and the current LRO-NAC derived albedo mosaic covers 
about 2% of the Lunar surface at 1 meter/pixel resolution. Com- 
pared to imagery obtained from more recent missions, the images 
captured on film by the Apollo missions have more noise artifacts 
and errors in camera position/orientation and exposure time, mak- 
ing the processing of this data significantly more challenging. 

The reconstructed albedo is shown in Figure 5 over the Clemen- 
tine mission image mosaic. The images captured by the Clemen- 
tine mission at 500 meters/pixel resolution were taken at very low 
incidence and emission angles (that is, both the Sun and the space- 
craft were rather high in the sky above the terrain) making these 
images have an average intensity close to the average albedo. Note 
that the transition between the reconstructed albedo and the back- 
ground Clementine data is practically seamless demonstrating that 
our approach can recover the average albedo intensity. However, 
the Clementine images do not describe a true albedo map since at 
each pixel the local emission and incidence angles are non-zero 
due to the irregular terrain shape. 

Our approach significantly removes shadow and illumination 
artifacts in the input imagery. Note the shadow removal effect in 
Figure 7 (top-left), and the contour of crater showing the ejecta 


Fig. 5. Reconstructed albedo of the Apollo zone on top of Clemen- 
tine imagery in Google Earth. 



Fig. 6. Error map for the reconstructed albedo. 

The Lunar albedo and the corresponding error map files can 
be accessed in GeoTiff tile format from 

http s : //by s s . arc . nas a. gov/albedo/albedo . html 


and 


https://byss.arc.nasa.gov/albedo/error.html 





Fig. 7. Image mosaic of the Apollo zone (left), and the recon- 
structed albedo (right). 



Fig. 8. Image mosaic of the Apollo zone (left), and the recon- 
structed albedo (right). 


respectively. The albedo in KML format (viewable in Google 
Earth) can be accessed at 

http s : //by s s . arc . nas a. gov/albedo/albedo . kml 

5. CONCLUSIONS AND FUTURE WORK 

This work introduces a method for photometric reconstruction of 
the Lunar surface from archival images taken during the Apollo 15, 
16 and 17 missions. The reconstructed area covers approximately 
16% of the Lunar surface at an unprecedented resolution of 10 me- 
ters/pixel for this coverage. The proposed method jointly estimates 
the surface albedo and camera exposure time, and refines the Lu- 
nar photometric parameters using a weighted least squares method. 
The reconstructed albedo and reconstruction error are made pub- 
licly available, and it is hoped that they will help planetary sci- 
entists and mission planners in current and future NASA missions. 
The methods introduced in this paper also enable future research in 


robust surface reconstruction and shape-from- shading techniques 
that will further increase the resolution of current terrain models. 
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